The criticality of the Hantavirus infected phase at Zuni 
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A preliminary analysis of the temporal evolution of a population of Peromyscus maniculatus 
infected with Hantavirus Sin Nombre is made. Ecological and epidemiological parameters are derived 
from the data, and they are used as inputs for the analytical model presented in [Abramson and 
Kenkre, Phys. Rev. E 66, 011912 (2002)]. A prediction of the critical carrying capacity and its 
associated critical mouse density is made, and the time series is analyzed under the light of these. 
It is found that the sporadic disappearances and reappearances of the infected phase correspond to 
the bifurcation predicted by the model. 
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This report contains a first attempt to test, in a quan- 
titative way, the predictions of the analytical model for 
the Hantavirus epizootic proposed by Abramson and 
Kenkre Q, 0- The analysis is based on ecological and 
epidemic data published by Yates et al. , consisting of 
a time series of mice populations at two sites near Zuni, 
New Mexico, covering monthly a period from 1995 to 
2001 (78 months). The population data are displayed 
in Fig. ^ where the total, susceptible and infected pop- 
ulations of Peromyscus maniculatus, the main host of 
Hantavirus Sin Nombre, are shown as N, N$ and Nj 
respectively. The period does not include the outbreak 
of 1993, when the virus was discovered to be the cause 
of the severe Hantavirus Pulmonary Syndrome (HPS). 
The great peak that can be seen around month 40 in 
Fig. n corresponds to a population explosion following 
the Nino event of 1997-1998, that was in turn followed 
by a small outbreak of AT/ and HPS in 1999 (see Fig-HJ. 
Observe also that the vertical scale is logarithmic to fa- 
cilitate the visualization of small values. It can be seen 
that Nj drops to several times, and the corresponding 
line is discontinuous at these values because of the log- 
arithm (not because of a lack of data). These are the 
sporadic disappearance of the infected phase, reported 
in d, 0, 0j an d that are accounted for as bifurcations 
controlled by the environmental carrying capacity in the 
model of Ref. p) . The assumptions of the present analy- 
sis are big but the results seem encouraging and support 
the analytical results of Ref. [l| . 

The basic model, mean-ficld-like and not extended in 
space, for the dynamics of the mice populations is: 

dNs = bN-cNs-^^-aNsN!, (1) 
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where b is the birth rate, c is the death rate, a is the 
contagion rate, and K is the carrying capacity. The in- 
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FIG. 1: Population of P. maniculatus at the Zuni capture site 
(New Mexico), from 3]. Total, susceptible and infected po_ 
ulations are shown, as indicated in the legend. Reference 
provides data for the total and the seropositive populations, 
here represented as N and Nz respectively. The susceptible 
population Ns is defined as N — Ni. 



terested reader may consult Refs. 0,0 for a detailed 
discussion on the motivation and the implications of this 
model, as well as for a discussion of infection waves ob- 
served in the spatially extended model. 

Let us make certain assumptions that will allow us to 
derive values for the parameters from the data. The first 
assumption is that 6, c and a are constant, determined 
by biological properties of the agents involved, and inde- 
pendent of anything that is changing in the system and 
that makes the populations grow and decline. This role 
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is reserved solely for K = Kit). This is, indeed, the same 
assumption made in Refs. (ll|2j, to analyze the bifurca- 
tion and to simulate a time dependent scenario. 

A further assumption is to use simplified models to 
extract the parameters from restricted regions of the 
graphs. For example, observe the best visible popula- 
tion explosion of N, from t = 30 to 40. It is a very well 



2 




FIG. 2: The connection between the environment, in par- 
ticular weather, and incidence of HPS becomes apparent 
in this figure, that shows an average sea surface tempera- 
ture (SST) of the East Pacific Ocean, average precipitation 
(PPT) in Southwest North America, and cases of HPS in 
the USA (total and in the Four Corners region, near Zuni). 
Sources: The SST has been calculated as an average of all the 
East Pacific buoys of the TAO (Tropical Atmosphere Ocean) 
Project (www.pmel.noaa.gov/tao); the PPT are the "South- 
west" values provided by the National Climatic Data Center 
(www.ncdc.noaa.gov); HPS cases are from Centers for Disease 
Control and Prevention (www.cdc.gov). 



defined exponential, and it strongly suggests the use of 



(3) 



that is, the linear part of the logistic model for N. So, 
the simplest way to calculate b — c is to suppose an ex- 
ponential growth in all the total population increases. 

If the processes of natural death and contagion can 
be supposed decoupled to some extent, we can extract a 
from the explosions of Ni (contagion is the only source 
of Ni) and c from its decays (since —cNj is the linear 
part of the sink of TV/). 

So, the third assumptions is that the decays of Nj can 
be supposed exponential: 



= -cNr 



(4) 



whence the parameter c can be obtained. 

Now, the only source of Nj is the infection, so in the 
explosions of iVj we suppose the simplified equation: 



d/V/ 

~~dT 



= aNiNs, 



(5) 



that constitutes the fourth assumption. 

The three Eqs. ©, 01 and © can be solved (this last 
only approximately) and the solutions fitted to the data. 



The exponential growth of iV is the nicest, because there 
are several relatively broad regimes of such behavior. The 
two other processes, decay and infection, are less clear 
cut, because both are derived from the same Ni(t), and 
also because the regimes are narrower. 
The results are the following: 



0.21, 



0.75, 



0.77, 



(6) 



in units of month -1 (c and b — c) and (mice x month) -1 
(a). It is imperative to take this values cum grano salts. 
Besides the assumptions mentioned above, there are fur- 
ther uncertainties. In the average a, for example, a value 
of 17 was discarded because it is so bigger than all the 
others that is inevitable to suspect some artifact of the 
time series that, at this point, cannot be clarified. Let us 
carry on, nevertheless, and from © and Ref. pj conclude 
that 



K,= 



a(b 



6 mice x month. 



(7) 



This value of K c is the main quantitative result of the 
present analysis. It is a prediction of the model, based 
on the numerical values of the parameters as can be esti- 
mated from the dynamics of the populations. How does 
it compare with features of the time series? If we observe 
that the equilibrium solution of the logistic equation sat- 
isfied by N, when K is constant, is N* — K(b — c), there 
is a "critical density" N c = K c (b — c) with which the 
actual density can be compared. When not growing or 
decaying, N(t) will be more or less at equilibrium during 
the periods that K(t) remains constant. N c = 2 mice, 
approximately corresponding to the calculated parame- 
ters, is shown as a black horizontal line in Fig. ^ 

Let us analyze the time series under the light of this 
result. At t < 15 the population is above critical, and 
there is, correspondingly, a positive infected phase. Dur- 
ing this regime, it is conceivable that K (t) has some time 
dependence, but that it keeps it above critical. Then, 
around t = 15, K (t) drops to some value below critical, 
and so does N*, and in consequence N(t) drops, trying to 
reach an equilibrium which is now below N c . The drop is 
not monotonous, there seems to be some discrete steps. 
What is the result? Shortly after the decline of N begins 
(the indication that K has gone subcritical) the infected 
phase begins to disappear sporadically. A few infected 
mice may be entering by migration, or marginal suscep- 
tible mice might be being infected, but it is clear that 
the infection is disappearing from the site. After t = 30 
a steady population explosion starts, indicating that K 
has increased, and the observation that the population 
grows beyond N c indicates that the system is supercrit- 
ical again. A recovering of the infected phase is to be 
expected. It takes time, however, just as in the model 
(see Fig. 2 in Ref.Q), and not before t — 40 do we see 
a positive Nj again. After this, the population remains 
above critical, so K(t) must be critical most of the time, 
and the infection persists. A brief excursion of N below 
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iV c at t w 70 might be the beginning of a new extinc- 
tion event, and indeed Ni reaches its lowest values since 
t = 40. But shortly after this the time series ends and 
the analysis can not be carried further. Observe that the 
drop in TV (so in K) takes place in 2001, a year that was 
particularly dry, as can be seen in the precipitation data 
in Fig. 

In summary, even if there are a number of uncertain- 
ties, the analysis fits nicely in the picture given by the 
model of Ref. p| , namely that there is a critical point and 
that the system lives close to it. It encourages to attempt 



a more detailed study, in particular involving bigger data 
sets of the same system. 
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